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Initial Kaluza-Klein fluctuations and inflationary gravitational waves in braneworld 

cosmology 

Tsutomu KobayashQ 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

We study the spectrum of gravitational waves generated from inflation in the Randall-Sundrum 
braneworld. Since the inflationary gravitational waves are of quantum-mechanical origin, the initial 
configuration of perturbations in the bulk includes Kaluza-Klein quantum fluctuations as well as 
fluctuations in the zero mode. We show, however, that the initial fluctuations in Kaluza-Klein modes 
have no significant effect on the late time spectrum, irrespective of the energy scale of inflation and 
the equation of state parameter in the post-inflationary stage. This is done numerically, using the 
Wronskian formulation. 
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I. INTRODUCTION 

The inflationary scenario predicts the existence of 
a gravitational wave background produced quantum- 
mechanically Q] , in much the same way as the mechanism 
of generating density perturbations. Detecting the infla- 
tionary gravitational wave background is a grand chal- 
lenge for LISA 0, BBO 0, and other missions 0]. It 
will open up a new window into the very early universe, 
and thus in principle the possibility of recently proposed 
braneworld models [j| can be tested using the inflation- 
ary gravitational waves. 

A number of studies have been done aiming at obtain- 
ing a clear picture of the generation and evolution of grav- 
itational waves in the Randall-Sundrum braneworld [aLl 

BHuamiiamQmiiapiiiiiaGaELEa- ^ 

Sitter inflation on the brane [6j is a special case where 
the perturbation equation is separable and exactly solv- 
able. It predicts a flat primordial spectrum, as in four- 
dimensional general relativity, but the amplitude is en- 
hanced at high energies. For more general inflation mod- 
els, the tilted primordial spectrum can be calculated from 
the corresponding four-dimensional one by a mapping 
formula Q . The late time evolution of gravitational 
wave perturbations at low energies is very close to that 
in conventional cosmology 0, liol recovering four- 
dimensional general relativity when the relevant length 
scales are lager than the bulk curvature scale £. 

What is much more interesting is the evolution of grav- 
itational waves during the high energy regime, £H ^> 1, 
where H is the Hubble parameter. In addition to the 
rather trivial effect of the unconventional background 
expansion rate, the nontrivial effect of mode mixing 
may not be negligible, converting a zero mode into a 
Kaluza-Klein (KK) mode and vice versa efficiently at 
high energies. To clarify the evolution of subhorizon 
modes during the high energy regime, several numer- 
ical studies have been done with different numerical 



schemes 0, 0, 0, 0, ^| . Most of them concentrated 
on the radiation-dominated phase after inflation, and 
very recently Hiramatsu 0] investigated consequences 
of a different equation of state parameter w. All of the 
above papers consider very naive initial conditions, ignor- 
ing initial abundance of KK fluctuations. However, since 
the gravitational wave perturbations are of quantum- 
mechanical origin, the most plausible initial configura- 
tion of the perturbations during inflation should include 
vacuum fluctuations in KK modes. Taking this point 
into account, in Refs. H,E1 Kobayashi and Tanaka have 
developed a numerical formulation using the Wronskian, 
by which one can discuss the quantum-mechanical gen- 
eration and subsequent evolution of gravitational waves. 
Following the previous works m the present paper 

we revisit the spectrum of inflationary gravitational wave 
background in the braneworld, focusing in particular on 
the effect of initial abundance of the KK fluctuations. 

The paper is organized as follows. In the next sec- 
tion we describe the background cosmological model, and 
then in Sec. Ill we review the Wronskian formulation to 
compute the spectrum of inflationary gravitational waves 
numerically. We present our numerical results in Sec. IV. 
Our conclusions are drawn in Sec. V. 



II. THE BACKGROUND MODEL 

We will work in the Randall-Sundrum- type 
braneworld p^ |. and the bulk metric is given by 
the five-dimensional anti-de Sitter spacctimc 



ds 2 = 



-dt 2 + 5ijdx l dx-i + dz 2 ) 



(1) 



where I is the bulk curvature scale. A cosmological so- 
lution can be described by a moving brane in the above 
static coordinates. The scale factor a(t) can be expressed 
in terms of the location of the brane z(t) as ait) = £/z(t), 
and it is governed by the modified Friedmann equa- 
tion [2! 
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where p is the matter energy density on the brane and 
(j = 6Afpj/£ 2 is the tension of the brane. The matter 
content is characterized by the equation of state param- 
eter 



W 



P 
P 



(3) 



Since the standard conservation law dp/ dr = —3H(p+p) 
holds on the brane, we have p oc a~ 3 ^ 1+wS> when w is 
constant. Using this and Eq. we obtain the scale 
factor as a function of the proper time r on the brane. 

The initial stage of the background model is assumed 
to be described by a de Sitter brane with a constant 
expansion rate Hi, and the final stage by a Minkowski 
brane. The two stages are connected smoothly by a 
Friedmann-Robertson- Walker (FRW) brane with a con- 
stant equation of state parameter w. We dub this phase 
as the "FRW phase" . The connection to the Minkowski 
phase is done at sufficiently low energies because we 
would like to focus on the possible high energy effects 
at the early stage of the FRW phase just after inflation. 
The construction here allows us to discuss the quantum- 
mechanical generation of gravitational waves during in- 
flation and their subsequent evolution during the FRW 
phase, and to see the final amplitude of the well-defined 
zero mode. In the previous paper ^7J orn y the radiation- 
dominated stage (w — 1/3) was considered, but in the 
present paper we are interested in other various values of 
w as well. 



III. FORMULATION 

A. Double null coordinates 

Let us consider gravitational wave (tensor-type) per- 
turbations. The perturbed metric can be written as 



ds 2 = — [-dt 2 + (5ij + ht^dx'dx 3 + dz 2 } , (4) 



where hij — hij(t,x, z) is the transverse-traceless metric 
perturbation. We decompose it into the spatial Fourier 
modes as 



(5) 



where M§ is the fundamental mass scale which is related 
to the four-dimensional Planck mass Mpi by (.(M^f = 
Mpy. The linearized Einstein equations give the Klein- 
Gordon- type equation for fa: 



\dt 2 



! +k 2 



d 2 3d 



ik = o, 



(0) 



dz 2 z dz 

and the junction conditions at the brane read 

™ a dafa\ hl . anc = 0, (7) 



where n a is the unit normal to the brane. From now on 
we suppress the subscript k. 

Following the previous papers 0, 0] we use the Wron- 
skian formulation to compute the spectrum of gravita- 
tional waves. Basic details of the formulation and numer- 
ical scheme are explained in Ref. @, but here we make 
a slight improvement. Since double null coordinates are 
convenient for numerical studies, our starting point is 



t — z, 



t 



(8) 



In this coordinate system, the trajectory of the brane can 
be specified arbitrarily by 



v — q(u). 



(9) 



To simplify this boundary trajectory, we make a coordi- 
nate transformation 

U = U, q{V) = v, (10) 

and then the position of the brane is given by 

U = V. (11) 

The following further coordinate transformation turns 
out to be useful for actual numerical calculations: 



X = - e -iJ u 

2 



T = --In (-—1 . 

2 V I 



(12) 
(13) 



Note that by an appropriate choice of the origin of the 
time coordinate t, we can always have U < and V < 0. 
The position of the brane is simply given by 



T = X. 



(14) 



We will perform numerical computations with equal grid 
spacing A in both T and X directions. The advantage of 
using these new coordinates is that it is easier to extend 
the computational domain far away from the brane in 
these coordinates than in the (U, V) coordinates which 
were used in the previous studies. The step size e in 
the U (V) coordinate increases with increasing \X\ (\T\), 
but we will adjust A so that the maximum of e is not 
so different from the step size used in Ref. ^j- Thanks 
to this choice of A, the (T, X) coordinates will not cause 
any problems in the resolution of the numerical analysis. 
Now the bulk metric is written as 



ds 2 = A 2 {T,X) \-2e- 2X/e q'(T)dTdX + 8 l:j dx l dx 3 



where 



A(T,X) := - 



2f 



q(T)+£e- 2X / e: 
q(T) := q{-le- 2T '% 



(15) 

(16) 
(17) 
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and a prime stands for a derivative with respect to the 
argument. The proper time r and the scale factor a on 
the brane are given respectively by 



The index v is related to the KK mass m as 



(18) 
(19) 



dr = a- y /2q'(T)e- T / i dT, 
a = A(T, T), 

and hence we have 

q'(T) = 2e- 2T / e (y/\ + PH* ~ £ff) 2 ■ 

Given the Hubble parameter H as a function of r, one 
can integrate Eqs. Ijl8(l and (121 )[l to obtain q as a function 
of T. 

In the (T, X) coordinates the Klein-Gordon- type equa- 
tion for <b reduces to 



(20) 



2d T d 



q(T) + £e- 2X / e 
+e- 2X ' l q'(T)k 2 ^ = 0, 
subject to the boundary condition 

[d T - d x ] 0| T=X = 



2e- 2X ' t d T - q'(T)d x 



B. Wronskian formulation 



(21) 
(22) 



Since the system we are considering has infinitely many 
degrees of freedom, consisting of a zero mode and a tower 
of KK modes, it is not advisable to set up the initial value 
problem for the evolution of each perturbation modes. 
Instead, we confine our attention to a single degree of 
freedom concerning the zero mode in the final Minkowski 
phase. The necessary information that determines the 
amplitude of the final zero mode can be picked up by 
using the Wronskian. Our Wronskian evaluated on a 
constant T hypersurface is defined by 



(0! • fa) := 2^ / dXA A {T,X) (fadxfa - 4>* 2 d x <t>i) , 

J — oo 

(23) 

which is independent of the choice of the hypersurface. 

Treating the graviton field <f> as an operator, in the 
initial de Sitter phase we expand it using the annihilation 
and creation operators, a n and fit, as 



$ = a cf)o + %4>l + dv (fi„0„ + a* 0*) 



(24) 



where <p n and its complex conjugate 0* are the mode 
functions in the de Sitter phase (n = for the zero mode 
and n = v for the KK modes). These mode functions 
form a complete orthonormal basis with respect to the 
Wronskian: 



h • 0o) = -(0o ' 0o) = !) 

^0 • <f>u) = (00 • 4>t) = °' 

i>n ' 0n') = 0, for n, n = 0, v. 



(25) 



m 2 = 



1 1 



(26) 



and v > 0. 

Similarly, in the final Minkowski phase we expand it 

as 

4> = A Va+Altp* + ^ dm [A mVm + Alcp*^, (27) 

where (p n and y>* (n = 0, m) are the mode functions in 
the Minkowski phase, and A n and are the annihilation 
and creation operators, respectively, of their correspond- 
ing modes. These mode functions also form a complete 
orthonormal basis with respect to the Wronskian: 

(<p ■ <po) = -(<Po ■ fl) = 1, (28) 

(<Pm ■ <Pm') = ~(<Pm ' <Pm') = S(m - m'), 
{ifO ■ <p m ) = (<Po ■ V*rn) = °; 

(ip n ■ ip* n ,) = 0, for n, n' = 0, m, 

The explicit form of the above mode functions is pre- 
sented in Appendix. 

Using the Wronskian, we now explain how to compute 
the spectrum of gravitational waves. We assume that ini- 
tially the gravitons are in the de Sitter invariant vacuum 
state annihilated by So and a„, 



a |0) = a„|0) = 0. 



(29) 



The expectation value of the squared amplitude of the 
zero mode in the final Minkowski stage is 



(i ^o + i^S) 2 |0>^^iV/ 



(30) 



where Nf :— (0| Aq Ao|0) is the number of created zero 
mode gravitons and we assumed that Nf 3> 1. In 
deriving Eq. I|30|) we used the commutation relation 

j4o,^4o =1 and the expression of the zero mode func- 
tion (|A1|I . The final power spectrum is then given by 



T(k) := 



47rfc 3 



(2tt)3 (M 5 ) 3 



—Nt = 



k 2 



TT 2 M 2 X 



N 



/• 



(31) 



Note here that by multiplying the factor 2/(Ms) 3 we ob- 
tained the squared amplitude of the metric fluctuation, 
rather than that of the canonically normalized field [see 
the definition (JSJ]. 

In fact, particle production is suppressed, Nf < 1, for 
k > ki := ciiHi where ai is the scale factor at the end of 
inflation. In the present paper we are interested in the 
modes that exits the horizon during inflation, so that the 
relevant wavenumber is smaller than fcj. 

The operator Aq can be projected out by making use 
of the Wronskian relations. Noting that the Wronskian 
is constant in time, we have 



Aq = ((f)- ip ) 



$) 



<P0)f = 

$) l a + / dv((j) u ■ $),a„ + h.c, (32) 



4 



where $ is a solution of the Klein-Gordon equation l|21() 
whose final configuration is the zero mode function ip 
in the Minkowski phase, and subscript / and i denote 
the quantities on the final and initial hypersurfaces, re- 
spectively. It is clear that final zero mode gravitons are 
created from the vacuum fluctuations both in the initial 
zero mode and in the KK modes: 



with 



(33) 



Correspondingly, the power spectrum (|31|l can be written 
as a sum of the two contributions: 



where 



Vkk 



V = V +Vkk, 



TT 2 M^ 



TT 2 M^ 



(34) 

(35) 
(36) 



Thus what we need to do is to solve the backward evolu- 
tion of the field <E> and to evaluate the Wronskian on the 
initial hypersurface in the de Sitter phase, which can be 
done numerically pjlH7|. 



IV. INFLATIONARY GRAVITATIONAL WAVES 
WITH INITIAL KK FLUCTUATIONS 

Let St be the primordial amplitude of gravitational 
waves from de Sitter inflation, which is given by |6j] 



Si 



2C 2 (£H l ) 



Hi 
2tt 



(37) 



where Hi is the Hubble parameter during inflation and 
the function C(£Hi) is related to the normalization of 
the zero mode (see Appendix). After the end of inflation 
the brane universe is dominated by a perfect fluid whose 
equation of state parameter is w, and the FRW phase is 
connected to the Minkowski phase when the energy scale 
of the universe becomes sufficiently low: H = Hq <C t^ 1 ■ 
For the modes that reenter the horizon during the low 
energy regime, namely, for the long wavelength modes 
with fc -C fc* where 



a*/£, 



(38) 



corrections to their evolution are very small, suppressed 
by e and £ 2 ln£ jl ITT| : the mode mixing effect is in- 
efficient and gravity on the brane is basically described 
by four-dimensional general relativity. Hence, for such 
modes the spectrum of gravitational waves is expected 
to have the same spectrum as the four-dimensional one: 



(fc < k < fc*), 



(39) 



n 4D 



l + 3w 



(40) 



where fco is the wavenumber associated with the horizon 
scale at the end of the FRW phase: fco '■— ^qHq. For 
the relatively short wavelength modes with fc > fc* , three 
non-trivial things should be considered: (i) the modi- 
fied background expansion rate [Eq. (J2J], (ii) the exci- 
tation of KK modes relevant at high energies, and (iii) 
the initial quantum fluctuations in KK modes Vkk- If 
the gravitational waves propagating in the bulk are ne- 
glected, namely, if the effects (ii) and (iii) are neglected, 
one would obtain an enhanced spectrum for k > fc*, 



fc 



with 



2 + 3w 



(41) 



(42) 



However, this evaluation will not be correct because of ef- 
ficient mode mixing at high energies. 1 The true spectrum 
will be somewhat different from the four-dimensional one 
and differed also from Eqs. (|41() and (35) . We would like 
to investigate how the balance of the above three things 
affects the spectrum. 

In Refs. [Hill 13 the effects of (i) the modified Fried- 
mann equation and (ii) the KK mode excitation have 
been addressed using the numerical formulations very 
different from ours. 2 The authors of Refs. [l3L fl5L Hq 
focused on the classical evolution of perturbations, while 
the initial conditions they adopt are naive, neglecting 
the initial abundance of KK fluctuations. In our Wron- 
skian formulation, the initial conditions are set quantum- 
mechanically and hence are the most plausible. It has 
been (incompletely) shown that the effect of the ini- 
tial quantum fluctuations in the KK modes is subdom- 
inant relative to the other two effects in the radiation- 
dominated phase (w = 1/3) ^?|]- I n that case, however, 
the dominant two effects cancel each other and conse- 
quently we have the same spectrum as in the conventional 
four-dimensional universe. 

In the previous paper |17| the authors failed in the 
precise evaluation of the fraction coming from the initial 
KK modes, 



TKK(k) 



Vkk 



Vkk 



(43) 



1 There are also braneworld models where one arrives at an inter- 
esting conclusion abo ut cosmological gravitational waves without 
mode mixing 23. 1261 . 

2 The formulation by Hiramatsu et at. \\'-t is ba sed on the Poincare 
coordinates, while Ichiki and Nakamura ll ■"! used a single null 
coordinate. The two results do not agree with each other for some 
unclear reason. However, t he W ronskian approach gives the 
result consistent with Ref. ll 'I . 
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FIG. 1: (color online). Contribution of initial Kaluza-Klein 
fluctuations. The equation of state parameter is w = 1/3 and 
the inflationary energy scale is given by IHi = 100 (lower 
line), tkk is plotted also for IHi = 10 (upper line). (The 
previous estimate in Ref. UtI should be compared with the 
lower line.) 
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FIG. 2: (color online). Spectrum of gravitational waves from 
inflation with IHi = 100, followed by a FRW phase with 
w = 0. The total power spectrum is shown by red circles, 
while blue crosses represent the contribution only from the 
initial zero mode. Dashed line indicates the four-dimensional 
result. 



due to the limited number of grids in the extra direction. 
Now using the new coordinates T and X, we are able to 
obtain much more precise values of tkk- Our refined re- 
sult is shown in Fig.^ This is different from the previous 
estimate by a factor of two or so. 

To make it clear how the spectrum behaves in a more 
general situation, let us move on to the models with the 
equation of state parameter other than 1/3. We have per- 
formed numerical calculations for three different cases: 
w = 0,2/3, and 1. The results are shown in Figs. I2H3I 
We find that in all the cases the modification to the expan- 
sion rate and the KK mode excitation are the dominant 
effects, and the initial abundance of KK fluctuations gives 
a subdominant contribution. In contrast to the radiation- 
dominated model, the two dominant effects do not cancel 
each other, yielding the spectrum different from the four- 



FIG. 3: (color online). Same as Fig.0 but w = 2/3. 
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FIG. 4: (color online). Same as Fig. |5J but w = 1. 



dimensional one for short wavelength modes. Assuming 
that the spectrum has the form 



V(xk n , k>ak if , 
where we take a ~ 3, we obtain 



(44) 



(45) 



From this result and the previous one (n ~ —2 for 
w = 1/3), we may deduce that the spectral tilt for short 
wavelength modes is given generally by 



2.6, 


(w 


-0) 


1.8, 


(w 


= 2/3) 


1.6, 


(w 


= !)• 



5 + 3w 
2 + 3w ' 



(46) 



This a gree s with the recent numerical result by Hira- 
matsu |16( because of the insensitivity of the final result 
to the initial KK mode contamination. For the modes 
with k < fc*, we can confirm that the spectrum is ap- 
proximately given by the four-dimensional one in all the 
three cases. 

Now we take a closer look at the contribution from the 
initial KK fluctuations. The fraction vkk is plotted for 
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FIG. 5: (color online). Fraction tkk for different inflationary 
energy scales, showing (from up to bottom) IHi = 10, 50, and 
100. The equation of state parameter is given by w — 0. 
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FIG. 6: (color online). Same as Fig. but w = 1. 

various values of the inflationary Hubble parameter Hi in 
Figs. El and (and also in Fig.^J. While the equation of 
state parameter is different (w — for Fig. El and w = 1 
for Fig. EJ), the plots share the same feature. For fixed 
k, tkk increases with decreasing £Hi, but, interestingly, 
the maximum of tkk is given by ~ 0.1 irrespective of the 
inflationary energy scale IHi. We find that this behavior 
can be summarized in the following universal relation: 

r KK (k) ~ r max f — J , (47) 

with 

W < 0{l), (48) 

p « 3, (49) 

and note that fcj is given by 

h = OiHi ~ (^)(2+3-)/(3+3») x K (50) 

Our numerical analysis indicates that the formula 147|) 
holds as long as the inflationary stage lies in the high 
energy regime. For low energy inflation (IHi <C 1) tkk is 



further suppressed, and we confirmed that by computing 
rxK for a model with £Hi = 0.2 and w — 1/3, showing 
that tkk ~ 10 -3 at k ~ ajilj. 



V. CONCLUSIONS 

In this paper, we have investigated the spectrum 
of gravitational waves generated from inflation on the 
brane. Our analysis generalize the previous work [TtI ] in 
that the present background cosmological model allows 
for the post-inflationary stage characterized by the equa- 
tion of state parameter w other than 1/3. We have paid a 
particular attention to the effect of the initial condition: 
the configuration of the perturbations during inflation in- 
cludes vacuum fluctuations in Kaluza-Klein modes. This 
was made possible by the use of the Wronskian formula- 
tion [1113, which is a distinctive point compared to the 
other related works. 

Our numerical analysis for the various inflationary en- 
ergy scale and equation of state parameter showed that 
the contribution of the initial KK fluctuations ?"kk:(&)> 
defined in Eq. (|4*3*|) . always behaves as Eq. <|T7Jl . implying 
that the initial KK fluctuations give rise to the subdom- 
inant effect on the spectrum. This is the main conclu- 
sion of the paper, bein g tr ue in the more general setup 
than the previous work [l7| • The present calculation cov- 
ers the limited range of the inflationary energy scale (up 
to IHi — 100), but we believe the results to be suffi- 
cient for concluding that also for higher energy scales, 
tkk is smaller than order of unity even at the small- 
est wavelength concerned (k ~ fcj). While the effect 
of the initial KK fluctuations can be safely neglected, 
the enhancement effect due to the modified Friedmann 
equation and the damping effect due to the excitation of 
KK modes are dominantly work, yielding the spectral tilt 
n ~ — (5+3u>) / (2+3ui) for short wavelength modes. This 
can be a observational signature of the braneworld model 
with nonstandard (i.e., w ^ 1/3) cosmological history. 

All the results in this paper are given numerically, 
but surely they will be great help for understanding 
the evolution of gravitational wave perturbations in the 
braneworld. Now the clues are in order: the spectral 
index is given by Eq. H46fl and it is insensitive to the con- 
tamination of the initial KK fluctuations, with their effect 
expressed in detail by Eqs. I|47|l - (|49|l . We will hopefully 
provide some analytic arguments complementary to the 
present numerical results in a future publication. 
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APPENDIX A: MODE FUNCTIONS 

In this appendix, we present an explicit form of 
the mode functions in the Minkowski and de Sitter 
braneworlds. 

In the Minkowski braneworld, the normalized zero 
mode function is given by 



1 



-ikt 



>2M 



while the normalized KK mode function is 

1 



with 

u m (z) := z 2 

and 



-iujt 



V2^~ Um[Zl 
mY\(rn€)3i(m,z) — Ji(m£)Y2(mz) 



to 



(Al) 



(A2) 



,(A3) 



(A4) 



In the de Sitter braneworld we introduce new coordi- 
nates (?7, £), which is related to (t, z) as 



t = ?7COsh£ + to, z — — T^sinh^, 



(A5) 



where to is an arbitrary constant. In (/?,£) frame the 
de Sitter brane is located at a fixed coordinate position 
£ = £t = constant, and the Hubble parameter on the 



brane is given by Hi = i 1 sinh£b- The normalized zero 
mode is 



Mv) = C(tHi) 

with 
C(x) := 



<2H \ V /•- ' ' 



-ikrj 



yjl + x 2 + x 2 In 



l + Vl + a; 2 



-1/2 



(A6) 



(A7) 



and the KK mode functions are found in the form of 

<t>v(v,0 = i>v{v)Xv{t), where 



: r 3/a e -^/2 ( _ f7) 3/2 jff W(_ fcf7 ) ) (A8) 



XAO = ^(sinhC) 2 PT^ (coshO 



C 2 Q_ 2 1/2+iv (coshO 



(A9) 



with 



Ci 



t{w) 



Y(b/2 + iv 

r(-iiz) 



TrC; 



r(«i/-3/2) 



c 2 = 



r(5/2-ii/) r(l + ii/) 

P -l/2+ W ( COSh ^) 



-1/2 



:i/2+ W ( COSh ^)' 



(A10) 
(All) 
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